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ABSTRACT 

We study the tensor perturbations in a class of non-local, purely gravi- 
tational models which naturally end inflation in a distinctive phase of oscil- 
lations with slight and short violations of the weak energy condition. We 
find the usual generic form for the tensor power spectrum. The presence of 
the oscillatory phase leads to an enhancement of gravitational waves with 
frequencies somewhat less than 10 10 ffe. 
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• Introduction: During the inflationary era infrared virtual gravitons 
are ripped out of the vacuum by the accelerated expansion of spacetime 
PQ. We have proposed a mechanism through which the back- reaction to 
this process can slow the inflationary expansion rate [2]. The idea is that the 
energy density of the newly produced gravitons comes not just from their bare 
kinetic energy but also from the interactions they have with other gravitons. 
Because gravity is attractive, these interactions act to slow inflation. If one 
imagines the process occurring on a closed spatial manifold such as T 3 it 
is easy to show that the energy density of this gravitational self-interaction 
grows as the exponential of twice the number of e-foldings, assuming all 
the gravitons on a 3-surface were in interaction with one another [3J. That 
would be make the manifold suffer gravitational collapse after only about 
14 e-foldings! Of course this estimate ignores causality; any given graviton 
only feels the potentials from the gravitons which are visible within its past 
light-cone. The actual growth of the interaction energy is therefore like the 
number of e-foldings, so the system requires something like 10 6 e-foldings of 
inflation to evolve itself to the point of gravitational collapse. However, it 
seems inevitable that this effect must eventually stop inflation if nothing else 
supervenes. 

The mechanism just sketched provides a wonderful way of resolving two 
perplexing problems: 

• Explaining why the dimensionless product of Newton's constant G and 
the cosmological constant A is about 10~ 122 ; and 

• Finding a natural model of inflation. 

Our solution to the first problem is to deny its premise: we do not believe A is 
actually small. It only seems to be small when one infers its value from the 
current expansion rate using classical gravity, which ignores the quantum 
gravitational screening mechanism. To be specific, we believe GA < 10~ 6 
is just small enough to justify using the semi-classical approximation, and 
that this is what started primordial inflation, without the need for a scalar 
inflaton. This gives a very attractive model of inflation because: 

• It dispenses with the need to assume that the initial configuration of 
the inflaton is homogeneous over more than a Hubble volume since the 
cosmological constant is homogeneous over all space; 
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• The long period of inflation required for back-reaction to build up dis- 
penses with the problem of having to fine tune the inflaton potential 
to make inflation last long enough; 

• The screening mechanism dispenses with the need to fine tune the 
inflaton potential so that inflation ends with nearly zero cosmological 
constant; 

• The fact that inflation is driven by the bare cosmological constant, 
rather than a scalar inflaton, dispenses with the need to invent a new, 
otherwise undetected field; and 

• The fact that this model contains only the single dimensionless param- 
eter GA means that it makes unique predictions, unlike scalar-driven 
inflation which can be tuned to give almost any result for the scalar 
power spectrum. 

However, there is a severe problem of tractability when it comes to using 
this model. The bare kinetic energy of inflationary gravitons is a one-loop 
effect, so the interaction energy we seek does not appear until two- loop or- 
der. Two-loop results are not easy to compute even for simple theories on 
flat space background; this one requires quantum gravity on de Sitter back- 
ground. In fact, the graviton 1-point function has been evaluated at two- loop 
order and it does seem to support the relaxation mechanism [I], but this 
computation required a year's labor and its interpretation is open to debate 
owing to the difficulty of formulating an invariant definition of relaxation [5]. 
Furthermore, at the same time as the two-loop effect finally becomes signif- 
icant, the effects of higher loops also become significant. So one requires a 
non-perturbative resummation technique to evolve into the late time regime 
at which interesting predictions can be made. 

We believe it may be possible to derive such a non-perturbative resumma- 
tion technique by extending the stochastic method which Starobinsky devised 
for the same purpose in scalar potential models [6j [7] . However, generalizing 
this technique to gravity is a difficult problem [8] . This paper is part of a par- 
allel effort [9] which is based on the idea of guessing the most cosmologically 
significant part of the effective field equations of quantum gravity. While 
there is no chance of guessing the full effective field equations, it might be 
possible to guess just enough to correctly describe the evolution of the scale 
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factor a(t) for a homogeneous and isotropic geometry, using what we know 
from perturbation theory about how the back-reaction effect scales. 

An important feature of any model which relaxes the cosmological con- 
stant is non-locality. Nonlocal models of cosmology have been much studied 
[3j [10] because they can avoid the problem that de Sitter must be a solution 
for any local, stable theory, and because non-local couplings between differ- 
ent times can ease fine tuning problems. In a previous paper [H] we proposed 
a phenomenological model which can provide evolution beyond perturbation 
theory. In one sentence, we constructed an effective conserved stress-energy 
tensor T^ u [g] which modifies the gravitational equations of motion: 

G^u = R^u - ^9^R = -kg^v + SirGT^g] . (1) 

and which, we hope, contains the most cosmologically significant part of the 
full effective quantum gravitational equations. 

What form to guess for T^ u [g] was motivated by what we seek to do, and 
by what we know from perturbation theory. We seek to describe cosmology, 
which implies homogeneous and isotropic geometries. When specialized to 
such a geometry the full effective stress tensor must take the perfect fluid 
form and we lose nothing by assuming that generally: 

Tp,\g] = {p + p)u ll u u + pg^ . (2) 

The relation between p[g], p[g] and u^[g] is heavily constrained by stress- 
energy conservation, but it is possible to specify one function for free. It 
turns out to be computationally simplest to take this free function to be 
the pressure [9j. We further require the pressure to be an ordinary function 
of some non-local scalar which grows like the number of e-foldings when 
specialized to de Sitter. The simplest choice for such a scalar is the inverse 
d'Alembertian acting on the Ricci scalar [3J. If the pressure is to grow the 
way we know it does from perturbation theory [UJ, and to eventually end 
inflation, then one is lead to the form [9]: 

p[g](x) = A 2 f[-GAX}(x) , X = i J2 , (3) 

1 Hellenic indices take on spacetime values while Latin indices take on space values. 
Our metric tensor g^ v has spacelike signature and our curvature tensor equals: R%u, v = 
r %, M + rQ w r % ~ (M <-> v). The initial Hubble parameter is 3H$ = A. We restrict 
our analysis to scales M = (A/87rG)i below the Planck mass Mpj = G~i so that the 
dimensionless coupling constant GA = ( M/Mpi ) 4 of the theory is small. 
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where the function / grows without bound and satisfies: 



f[—GAX] = -GAX + 0[{GA) 2 \ , (4) 
and where the scalar d'Alembertian: 

' dJg^V=gd v ) , (5) 



is defined with retarded boundary conditions. The induced energy density 
p[g] and 4- velocity u^[g] are determined, up to their initial value data, from 
stress-energy conservation: 

D^T^ = . (6) 

The 4-velocity was chosen to be timelike and normalized: 

g^u^ = -1 u»uw = . (7) 

The homogeneous and isotropic evolution H of this model - using a com- 
bination of numerical and analytical methods - revealed the following basic 
features: El 

- After the onset and during the era of inflation, the source X(t) grows while 
the curvature scalar R(t) and Hubble parameter H(t) decrease. 

- Inflationary evolution dominates roughly until we reach a critical point X cr 
defined by: 

1 — 87rGA f[—GAX cr ] = . (8) 

- The epoch of inflation ends close to but before the universe evolves to the 
critical time. This is most directly seen from the deceleration parameter since 
initially q(t = 0) = — 1 while at criticality q(t = t cr ) = +|. 

- Oscillations in R(t) become significant as we approach the end of inflation; 
they are centered around R = 0, their frequency equals: 



u = GAH J72nf' cr , (9) 



2 The line element in co-moving coordinates is ds 2 = —dt 2 + a 2 (t) dx-dx. In terms of the 
scale factor a, the Hubble parameter equals H(t) = aa~ l and the deceleration parameter 
equals q(t) = -add^ 2 = -1 - H H~ 2 = -1 + e(t). 

3 In [9] , our analytical results were obtained for any function / satisfying (f4]) and growing 
without bound, our numerical results for the choice: f(x) — exp(x) — 1. 
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and their envelope is linearly falling with time. 

- During the oscillations era, although there is net expansion, the oscillations 
of H(t) take it to small negative values for short time intervals - a feature 
conducive to rapid reheating; those of H(t) take it to positive values for 
about half the time; and, those of a(t) are centered around a linear increase 
with time. 

Furthermore, this simple model can be improved to account for the re- 
quirement to naturally produce negative pressure of the correct magnitude 
during the current epoch |12j . However, for the purposes of this paper such 
an improvement is inconsequential since we shall be concerned with the pe- 
riod about the end of inflation where the simple model and its improvement 
are indistinguishable. 

A novel feature of this class of models is the existence of an oscillatory 
regime of short duration which commences towards the very end of the infla- 
tionary era. During this period H(t) is positive about half the time, which 
represents a violation of the weak energy condition. Such a violation can- 
not occur in classically stable theories [13] but it can be driven by quantum 
effects of the type we seek to model without endangering stability. 

A simple, fully worked-out example is provided by a massless, minimally 
coupled scalar with a 4 potential on non-dynamical de Sitter background. 
For that model, inflationary particle production tends to push the scalar 
up its potential, which of course increases the vacuum energy and leads to 
a violation of the weak energy condition. That has been confirmed in a 
fully renormalized computation of the expectation value of the model's stress 
tensor at two-loop order [13]. That there is no instability was confirmed by 
a fully renormalized computation of the scalar self-mass-squared at two-loop 
order, which was then used to solve the linearized effective field equations 
[TB] . As might be expected, pushing the scalar up its potential makes it 
develop a small mass, which actually makes the system more stable, not 
less. And a fully non-perturbative analysis by Starobinsky and Yokoyama 
confirms that the system approaches a static limit [IB] . 

Our own model is an attempt to model the most cosmologically significant 
features of the inflationary production of gravitons, so it shares many features 
with the simpler 4 model. Of course it is quite a bit more difficult derive 
comparably powerful results for quantum gravity. However, we have been 
able to show that there are no tachyonic modes [T7] . 

The oscillatory phase and the associated weak energy condition viola- 



5 



tion is a very distinctive feature of our model and cannot occur in classical, 
scalar-driven inflation. The purpose of this study is to determine whether 
this oscillatory regime leaves its signature on the observable tensor power 
spectrum. @ We shall, therefore, obtain the amplitude and frequency of two 
kinds of gravitational waves and examine their evolution under the expansion 
history that this class of models predicts. The first kind of waves is now of 
cosmological scale and originated during inflation while the second kind was 
on the verge of experiencing first horizon crossing when the epoch of oscilla- 
tions began. Because we shall be relating scales from the very early universe 
to current measurements, we first focus on presenting the basic equations 
and relevant relations, and then we apply them for our purpose. 

• The Set-up: The analysis of tensor perturbations in this class of 
models is much simpler than that of scalar perturbations |17J . The reason 
is that - unlike the case of scalar perturbations - the non-local nature of 
the model does not alter the basic equation which tensor perturbations hj^ 
satisfy at linearized order: 



d 2 d V 2 



dt 2 w dt a 2 (t) 



K (t, x) = o . (io) 



Therefore, up to sub-dominant corrections coming from the exact form of the 
mode functions before and after first horizon crossing, the resulting power 
spectrum A\ will have the usual form: 

1 fi 

A 2 (A;) ~ -GH\t k ) , (11) 

7T 

where the Hubble parameter H is evaluated at the time t k of first horizon 
crossing of the mode with wavenumber k: 

k = H(t k )a(t k ) . (12) 

Moreover, the tensor spectral index n T is defined as: 



4 The more complicated analysis of the scalar power spectrum that this class of non-local 
cosmological models predicts has been done elsewhere [17] . 
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and for the power spectrum (TTTi) equals: 

2e{t k ) 



1 - e(t h 



(14) 



where the last relation assumes that e(tfc) <C 1. It is apparent that knowl- 
edge of the relevant scale factor a(t) suffices to compute the tensor power 
spectrum and spectral index. 

• Assumptions about the Expansion History: For the purposes of 
this paper we divide cosmological history into three epochs: 

- Primordial Inflation. The most convenient time parameter for the epoch 
of primordial inflation is the number N of e-foldings before criticality: 

a{t) = a CT e~ N . (15) 

The important cosmological parameters during this phase are |17j : 

H\t) ~ i w 2 (4Af + |) , (16) 

*> - ' (17) 

At the end of inflation the scale factor is about a cr and the Hubble parameter 
is about uj. 

- Oscillations. The distinctive feature of our model is the epoch of oscilla- 
tions. The most convenient time parameter during this era is the co-moving 
time after criticality: 

At = t-t cr . (18) 
The important cosmological parameters during this phase are [9]: 



a(t) ~ a cv luAt + 1 + \/2 (cos(u> At) - 1 



H(t) 



UJ 



1 - y/2 sin(o; At) 



UJ 



At + (1 - V2) + V2 cos(a; At) 



(19) 
(20) 



V^[wAt+(l- V2) 1 cosfo; At) + 3 - 2^2 sm( u At) 
<t) * —± -U— '—— 2 • (21) 



1 - y/2 sin(uj At) 
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The epoch of oscillations is terminated by the flow of energy density to the 
matter sector from the vast reservoir of super-horizon scalar modes, all of 
which begin to oscillate with frequency u. As explained in our analysis of 
scalar perturbations [IT] , we believe these oscillations should lead to very 
rapid reheating. Let us call the number of oscillatory e-foldings AiV. Then, 
at the end of the oscillations era: 

a ~ a cr e AN , H ~ we" M . (22) 

Before proceeding, it is worth reviewing the argument - in section 6 of 
[TT] for why the oscillatory phase should lead to rapid reheating. First, 
note that causality imposes absolutely no obstacle to the decay, eventually 
into radiation, of a super-horizon mode which is oscillating at a frequency 
comparable to the Hubble scale. Indeed, the usual mechanism of reheating 
relies on precisely such a decay of the inflaton zero mode, which has infinite 
frequency. @ Second, note that the quantum gravitational back-reaction to 
inflation, which our non-local field equations seek to model, requires a very 
long period of nearly de Sitter inflation before enough inflationary gravitons 
have been created that their accumulated interaction is enough to slow in- 
flation. One can estimate the number of inflationary e-foldings to be about 
(GA) -1 > 10 6 . Third, this long period of inflation means that the number 
density of modes which have undergone first horizon crossing is staggering 

«~^exp(A) > # 3 xl0 106 . (23) 



3vr 2 ^GA' 

Finally, recall that in our model, all of these super-horizon modes begin os- 
cillating at about the same time. In the face of a number density such as 
(j23|) it makes little sense to attemp to estimate the rate at which the distur- 
bance of a single oscillating mode would lead to the production of relativistic 
matter. As long as that rate is non-zero - and the universal character of the 
gravitational coupling ensures it is non-zero - then the vast number of modes 
which participate must make reheating almost instantaneous. 

- ACDM. The ACDM cosmology after the epoch of oscillations is stan- 
dard, and we do not require explicit forms for the three geometrical param- 
eters. To compare quantities from the first two eras with their redshifted 
descendants at present time it is useful to express the energy density pr at 



'See section 5.5 of the text by Mukhanov [18] . 
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the onset of the ACDM epoch in terms of the reheating temperature Tr and 
the number n ~ 10 3 of relativistic species: 

3cV e-^« it 2 (k B T R )> 

SnG "" X 30lfeF- (24) 

The current energy density p now can be written in terms of its tiny radiation 
fraction fi r » 8.5 x 1CT 4 and the corresponding temperature T now m 2.726-K" 
of that radiation 



3c 2 H now 2 7r 2 (^ijT now ) 4 (ot^\ 

Pnow " ~^g~ ~ n~ r x 30 (hcf • 1 j 

Dividing (|25l) by (124)) gives a relation between current conditions and those 
prevailing at the end of inflation: 

-j^) {—) • (26) 

We define N now as the number of e-foldings from criticality to the present 
time. Using the relation ( 12 6 p and: 



T n n e AN 

- 1 now U R u cr c 



Tr 

we see that N now equals: 



(27) 



N now ~ AN + In 



T, 



R 



T 



2 



CO 



+ -AN - -InfenfiJ . (28) 
2 4 L r J 



We shall later show that the measured value of the scalar power spectrum 
A^ , and the current limit on the tensor-to-scalar ratio r, together imply the 
restriction u < 10 55 if now . Hence we conclude: 



iVnow < 63 + - AN . (29) 

The argument that u < 10 55 if now results from comparing expression ( fill 
with the current bound on the tensor contribution to the quadrupole moment. 
To make this comparison we must solve the following problem concerning the 
relation between late times and early times: 

01. Given a physical wave number K now at the current time, find the e- folding 
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iVhor when it experienced first horizon crossing during inflation. To solve this 
problem, we first use the horizon crossing condition (IT21) to express K now in 
terms of N^ or : 



k 



a(t k ] 



a(t k ) a c 
x — — - x — 



-uj Jan + - x e~ Nhor 

3 V 3 



x e 



(30) 
(31) 



Now invert (l3Tj) to solve for A^ or : 

1 



hor 



In 
1 



In 



- AT, 



now 



now "I - 2 ^ 



4 ,r 4 

-AT H 

9 27 



c 2 AT 2 



2n fL 



+ 



9 



ln( 



c 2 K 2 



(32) 
(33) 



For the £-th partial wave contribution to the anisotropies of the cosmic ray 
microwave background, the corresponding number A^ of e-foldings before the 
end of inflation is: 



A'now ~ 2 ^ 



(34) 



In ( -) + - In ( -^—) - -AN + - In (in n r ) + - In -In ( 



l 



The restriction uj < 10 55 H now then implies: 

N f < 65 



-ln(-) - -AN . 

2 V2/ 2 



(35) 



We cannot hope to detect a signal outside the range 2 < t < 100, so the 
interesting values of A^ lie within a band of only four e-foldings. 

We now deduce the restriction on oj coming from the measured value of 
the scalar power spectrum A\ [19J: 



A\{h 



2.44 x 10 



-9 



k = 0.0002 (MpcY 1 



(36) 



6 The inversion was done under the assumption that the Hubble parameter at the end 
of inflation is much bigger than its present value: ui H^J^ 1. 



10 



and the 95% confidence bound on the tensor-to-scalar ratio: r(k ) < 0.22 
j. Employing expressions (ITT]) and (I2"0"j) we get: 

Al(k ) = r(k )A 2 n (k ) * ^Gc 2 (4iVo + ^) (37) 

y7T v o ' 



< 



0.22 



2.44 x 10" 9 . (38) 



Now the wave number ko and its associated number of e-foldings Nq cor- 
respond to the I = 2 partial wave so that - under the asumption that the 
Hubble parameter at the end of inflation is much bigger than its present 
value (uH^ w ^> 1) and that the duration of the oscillations era is very short 
(AN < 10) - equation (I3~4|) implies: iVo = N i=2 > 60. Thus, expressions 
( I37ll38|) reduce to uy/G <2x 10~ 6 and when we convert to Hz we get: 



co < 2 x 10" 6 ^ « 3.7 x 10 37 Hz =► u < 10 55 H QOW , (39) 

where we used H now « 3.2 x 10 _18 ife for the current value of the Hubble 
parameter. 

At this stage it is natural to consider a second problem which is in some 
ways the inverse of the first: 

02. Given a physical wave number K N from the epoch of inflation, find its 
physical wave number now. To achieve this, we express the current physical 
wave number in terms of Kn and N: 

^now = — = X — X — = K " X ^ X 

Sow a\tj O cr C^now 



K%H novl (nVtr^i ^.tv-IaAt ^q-j 



where in the last step we used ( 128]) for iV now . An important special case is 
the oscillation frequency / pea k for the wave vector — ^ at N = . In 
this situation, f l4"U|) gives: 



!^^U— ) « 2^5= (41) 
2tt \a now J Ko= ± 2vr 

Imposing the restriction uj < 10 55 if now and using the current value of H now « 
71 kms~ 1 (Mpc)- 1 « 3.2 x 10" 18 ife JTS], implies: 

/peak < (l0 9 /fe)e-5 A ^ . (42) 
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Therefore, the late time descendant of the gravitational waves produced by 
the phase of oscillations are unobservable in the cosmic microwave back- 
ground and have relevance only for high frequency direct detectors. 



• Overview of Gravitational Waves in the Oscillating Regime: 

In terms of the mode functions u(t, k) the basic equation (fTUl) takes the form: 



k 2 

u + 3H u + —u 

or 



. 



(43) 



We do not possess exact forms for the two, linearly independent solutions 
during the oscillatory regime. Even if we had these solutions, we would not 
know the linear combination of them that gives "the" mode function u(t, k), 
which we define to be the coefficient of the annihilation operator in the free 
field expansion of the graviton. It makes sense to first develop a reason- 
able approximation for the linearly independent solutions and then consider 
which combination of them occurs in the actual mode function u(t,k). In 
approximating the solutions it also makes sense to first include the effect 
of the overall linear expansion - for which exact solutions exist - and then 
numerically superimpose the effect of the oscillations. 



• The Case of Linear Expansion: During the oscillatory epoch the 
scale factor f[T5j) consists of a linear expansion plus an oscillatory term which 
causes the Hubble parameter f l2"U|) to become negative for brief periods. Be- 
cause we wish to quantify the potential enhancement from these periods of 
negative H(t), it is useful to factor out the behaviour that would arise from 
the linear growth, without the oscillatory term: 



a(t) 



1 + uAt 



(44) 



Then, the Hubble parameter H can be expressed in terms of the scale factor 
a as follows: 

H(t) = = ^ • (45) 

w 1 + uAt a(t) y ' 

The canonically normalized, Bunch-Davies mode function for the linear ex- 
pansion is [20J: 



u(t, k) 



1 



2 Jc 2 k 2 



u 2 a 2 



a(t) 



exp 



\ 



l k 2 



ui 2 a 2 



1 In 



a(t) 



• (46) 
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Since the product H(t) x a(t) is constant, there is no horizon crossing during 
linear expansion. Modes which are sub-horizon at criticality (k > u>a cr ) 
remain sub-horizon, and modes which are super-horizon at criticality (k < 
ua cv ) also remain super-horizon. One can see from expression (j46p that sub- 
horizon mode functions oscillate and redshift, with the period of oscillation 
also redshifting. By contrast, super-horizon mode functions fall off like: 



u 



super 



(t,k) 



a(t) 



i±Ji 



(47) 



• Initial Conditions: We still have to include the effect of oscillations, 
which can only be done numerically. That defines a mode function u(t, k) 
which obeys equation (f4"3l with the full oscillating geometry fll9H20p . We 
construct these mode functions to agree initially with those of the fictitious 
phase of linear expansion: 



u(t cr , k) — u(t cr , k) 



u(t cl[ , k) — u(t CT , k) . 



(48) 



The actual mode function u(t, k) - by which we mean the coefficient of the an- 
nihilation operator in the free field expansion - is neither u(t, k) nor u*(t, k), 
but rather a linear combination of the two solutions: 



u(t,k) = au(t,k) + f3u*(t,k) . 



(49) 



We can solve for the combination coefficients in terms of the values of u(t, k) 
and its first derivative at criticality: 



a 







-i a„ 



-i a" 



u(t cr , k) u (t cr , k) ii(t CT , k) u (t cr , k) 
u(t CT , k) u(t cr , k) - u(t cr , k) u(t cr , k) 



(50) 
(51) 



Although we do not know precisely what these values are, some reason- 
able guesses can be made. For example, a far super-horizon mode, which 
experienced first horizon crossing A^ or e-foldings before criticality, should 
have: 



u(t CT , k) 
ii(t cr , k) 



H 



/2F 

U N hol 



k 



2 r 



W \H Nhm a c 



1 + 



ik 



H Nhoi a CT _ 



(52) 
(53) 



7 The Wronskian of the two solutions is: 
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t 



2 4 6 8 10 12 14 



Figure 1: The enhancement factor Q versus co-moving time (in units of u> x ) for a 
super-horizon mode with ck = 0.01 x wo cr . 



Recall that the potentially observable modes in the cosmic microwave back- 
ground correspond to A/h or ~ 60, which implies: 

k - lO" 26 (54) 



H Nhoi a. 



There is absolutely no point in retaining such small numbers. So during 
the oscillatory phase after criticality, the mode function of a cosmologically 
observable wave number would be unchanged from ( |52|) . for all practical 
purposes. That had to be true because, for far super-horizon wave numbers, 
(l4"5|) simplifies to: 

u(t,k) + 3H(t)u(t,k) « . (55) 

and u(t, k) = constant remains a solution - independent of a(t) - for as long 
as it is valid to neglect the last term of ( )43l) . 



• Gravitational Waves Enhancement: Now let us consider the ef- 
fect of adding the oscillatory term to the scale factor. It is intuitively obvious 
that one gets a significant response at resonance; then, the natural time scale 
of the mode function is close to the inverse of the oscillatory frequency u. 
Whether or not this occurs depends upon two things: the wave number k 
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and the values of u(t, k) and ii(t, k) at the start of the oscillatory period. The 
reason the initial condition matters is that there are always two, linearly in- 
dependent solutions to the mode equation and they can have vastly different 
natural time scales. There are three interesting wave number regimes: 

* The far super-horizon, with ck <C ua cr . 
From expression (|47p it is evident that, without the oscillatory term, super- 
horizon modes fall off with time scales: 



For ck <C uia cr one of these is - within a factor of two - close to u while the 
other is vastly longer. Numerical analysis shows - see Figure [TJ - that the 
oscillations amplify the solution with the shorter time scale by about a factor 
of four. As might be expected, the solution with the longer time scale ex- 
periences no significant amplification. Because the natural initial conditions 
( 1521153]) imply the mode enters the oscillatory epoch almost entirely in the 
long time scale solution, the effect is that far super-horizon modes experience 
no significant enhancement from the oscillation. 

* The far sub-horizon, with ck ^> ua cr . 

Far sub-horizon modes also receive no substantial enhancement, but for a 
different reason. For far sub-horizon modes the natural frequencies of both 
solutions are about ck a" 1 which is much bigger than u, so neither solution 
experiences much enhancement and it does not matter much what the initial 
condition is. 

* The near-horizon, with ck « ua CT . 

As one might expect, it is the near horizon modes which experience the 
greatest enhancement. Figures [2] and [3] present numerical simulations for the 
case of ck = ^ua cr , giving the ratios of the actual mode functions - evolved 
with the oscillatory term - compared with the solution ( 146|) which starts from 
the same initial condition but is evolved without the oscillatory term. | In 
the near horizon regime one expects both solutions to be present with about 
the same amplitude, so a reasonable estimate of the total enhancement is 

8 The enhancement factor Ql is associated with the real part: Ql = Re[u(t, k)] 
Re[u(i, k)] while Ql with the imaginary part: Q2 = Im[u(t, k)] Im[w(t, fc)]. 



-1 




(56) 
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Figure 3: The enhancement factor Q2 versus co- moving time (in units of to x ) for a 
near-horizon mode with ck = 1.1 x usa CI . 
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by adding the two solutions in quadrature and taking the ratio with, and 
without the oscillatory term: 

\u(t, k) | 

Q = • (57) 

From Figure H] one can see that the enhancement factor is about Q ~ 10. 
By comparison with the Hubble parameter - see Figure [5]- we see that almost 
all the enhancement derives from the first oscillation. This is important 
because, as we explained, the phase of oscillations is likely to be short. The 
fact that almost all the enhancement occurs during the first oscillation means 
that the effect is reliable even if only one oscillation occurs. 

Figures [6] and [7] give the total enhancement factor for the cases of ck = 
2ua cr and ck = 5ua CT , respectively. Note that the enhancement falls off 
rapidly as one moves away from resonance. Note also that the enhancement 
factor oscillates. 

Modes which are slightly super-horizon are quite similar to those which 
are slightly sub-horizon. Figure |8] gives the total enhancement factor for 
the case of ck = |wa cr . However, decreasing the wave number much more 
rapidly reaches the factor of four enhancement which is concentrated on the 
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H 




Figure 5: The Hubble paramemter H versus co-moving time (in units of w). The first 
period of H < coincides with the largest growth of Q in Figure [4] . 

solution that is not likely to be present. Figures I§1 and [TU1 give the behaviours 
for ck = \uja cr and ck = \uja cr , respectively. 

• Enhanced Waves Energy Density and Frequency: It remains 
to estimate the current energy density and frequency of gravitons which are 
produced during the epoch of oscillations. Suppose we regard the enhance- 



ment factor as Q = 10 for modes within the range 



< ck < 



i;tUa r 



and zero outside this band. This is superimposed on the mode functions ( |46|) 
of linear expansion. A reasonable estimate for the extra physical energy - 
above the 0-point value of ^hck x a _1 (t) - in a single wave vector k within 
the band of enhancement is: 



E(t CT , k) 



h&k 



.2,2 



Q u(t CT , k) 



Q 2 hc 2 k 2 



2ualJ\l - f-£*- 



(58) 
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Figure 6: The enhancement factor Q versus co-moving time (in units of for a 

near- horizon mode with ck = 2 x u;a cr . 
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Figure 9: The enhancement factor Q versus co-moving time (in units of co x ) for a 
near-horizon mode with ck = 0.5 x ua a . 
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Figure 10: The enhancement factor Q versus co-moving time (in units of oj 1 ) for a 
near-horizon mode with ck = 0.2 x uo cr . 



The associated energy density comes from integrating the extra physical 
energy (158]) over all the modes whose wave numbers are within the band of 
enhancement: 

Pgw(^cr) ~ 4f J 7^jT^ 0(^ua CI - ck) 9(ck - ^ua c ^j E(t cr ,k) (59) 
~ : 9 • (60) 

47r 2 c 3 V / 

Perhaps of more relevance for gravity waves detectors is the amount of energy 
density per unit frequency since such detectors are only sensitive in certain 
frequency bands. We can estimate the energy density per angular frequency 
at the critical time by simply not performing the radial integration over k in 
expression ( l60l) : 

dp gw (t cr ) 100 hw* (g) 2 (^) 4 
dck 47r 2 c 3 a rr . /h _ (shIW\ 



war- 



To convert this to the current epoch we note that these gravitons are sub- 
horizon after the phase of oscillations, so their energy density redshifts like 
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radiation: 



*(^now) ^ ^ Pgw(^cr) 



(62) 



We also take note of the relations between the ordinary (not angular) current 
frequency and peak frequency: 



fn 



ck 



27ra r , 



uja r 



peak 



27ra r 



(63) 



and we divide out a factor of the current critical density to obtain the fraction 
in gravity waves: 

— Pgwl^now/ • 



8ttG 



(64) 



The final result is: 
dVL 



"gw 



./no 



/peak 



/200\ /hGu 2 \ / lo \2 / a cr \' 

( — )( —)(jT- ) (~ ) 

^ II nnw u now 



\2 (fn 



X 



10' 



v /peak ' 



' /now \2 I 

. f , / I 
J peak 



(65) 



Using the restrictions ( 1291 1391) gives 



gw 



d( 



./no 



/peak ' 



(4 x 10" 



10 



Q_\2 ( /now \4 

1 - (^L) 2 | 

/peak 



(66) 



A rough model for the Q-factor supported by our numerical analysis is: 



Q 

— ~ exp 
10 P 



\ /peak 



(67) 



and it is for this choice of the Q-factor that Figure [TT] gives a plot of 
[dtt gw /d(j^)] as a function of [/now //peak]- 

From Figure [11] it is evident that the signal is highly peaked at the fre- 
quency /peak ~ 10 9 ffe, and is negligible at significantly different frequencies. 
It would be challenging to detect gravitational radiation at such high fre- 
quencies but detectors in that range have been proposed [21]. As noted in 
the text, the phase of oscillations does not affect modes which experienced 
first horizon crossing more than a few e-foldings before the end of inflation. 
The wavelength of our effect is A pea k = c/"^ > 0.3m, whereas the smallest 
scale feature which is currently observed in the cosmic microwave radiation 
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0.5 1.0 1.5 2.0 /peak 

Figure 1 1 : Fraction of the current energy density in gravity waves from our signal per 
frequency. 

is about 10 22 m [22J! Our model does not change either how matter couples 
to gravity or the propagation of linearized gravitons, so it has no effect on 
the spin-down rate of the binary pulsars. The gravity waves we predict will 
certainly distort how pulsar light propagates, but the short wavelength again 
seems to preclude a detectable effect. LIGO is not sensitive above frequencies 
of 7000Hz, which is far too low. The situation is even worse with LISA's 
high frequency cutoff of 0.1Hz. 
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